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Abstract 

We find approximate solutions to the renormalization group equation which gov- 
erns the quantum evolution of the effective theory for the Color Glass Condensate. 
This is a functional Fokker-Planck equation which generates in particular the non- 
linear evolution equations previously derived by Balitsky and Kovchegov within 
perturbative QCD. In the limit where the transverse momentum of the external 
probe is large compared to the saturation momentum, our approximations yield the 
Gaussian Ansatz for the effective action of the McLerran-Venugopalan model. In 
the opposite limit, of a small external momentum, we find that the effective theory 
is governed by a scale-invariant universal action which has the correct properties to 
describe gluon saturation. 
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There has been significant progress in recent years towards understanding the hadron 
structure in the high-density regime at small x, which is the regime relevant for high- 
energy scattering. As x decreases, the gluon density in the hadron wavefunction grows 
faster than the quark density [|l], giving rise to a multiparticle gluonic state with high 
occupation numbers. In this high density environment, important non-linear phenomena 
are expected, which should eventually lead to saturation [0, ^ ^, ^ 0, that is, to a limi- 
tation of the maximum gluon density per unit of phase-space. Linear evolution equations 
like BFKL [|l| do not take into account the rescattering among the produced gluons, and 
predict an exponential growth of the gluon distribution with ln(l/x). This leads to cross 
sections which in the high-energy limit violate the Froissart unitarity bound. If parton 
distributions saturate, then there is a natural resolution to this unitarity problem. 

Non-linear generalizations of the BFKL equation have been recently proposed within 
different approaches, with results which appear to be consistent with each other H, [T^. 



By using the operator expansion for high-energy scattering in QCD, Balitsky has derived 
a set of coupled equations for the evolution of Wilson-line operators p]. This is an 
infinite hierarchy of coupled equations for n-point operators, which however decouple in 
the large limit, with N the number of colors. In this limit, the equation satisfied by 
the 2-point function has been independently derived by Kovchegov within the dipole 



model of Mueller [O]. Solutions to Kovchegov's equation have been investigated in Refs. 



1^, |T2|, ITsf. In a recent paper [Q, Weigert has shown that Balitsky's equations can 



be summarized as a functional evolution equation for the generating functional of the 
Wilson line operators. Together with Leonidov, we have subsequently shown that 



this is precisely the renormalization group equation (RGE) which governs the quantum 



evolution of the effective theory for the Color Glass Condensate (CGC) |jT5|, [16 



In the description of saturation, a natural intrinsic momentum scale appears, the 
saturation momentum Qg ^, |T^, [1^, ^ . This is the scale below which the non-linear 



effects are expected to slow down and eventually saturate the increase of the gluon density 
with l//cl, where k± is the gluon transverse momentum^. In this paper, we shall find 
approximate solutions to the RGE for the CGC in two limits: at large momenta k± 3> Qs, 
and at small momenta fc_L -C Q^, where k± is the momentum scale at which we compute 
correlation functions, (a) In the first case, our results reproduce the Gaussian Ansatz 
for the "effective action" (strictly speaking, a probability distribution for a classical color 
source; see below) which has been postulated in the McLerran-Venugopalan (MV) model 
1^. We shall find improvements over the previous discussions of this model [^, [l^, |l8l in 
the sense of including the quantum evolution of the width of the Gaussian and properly 
treating the transverse non-locality. We shall also clarify the limits of validity of this 
model, and show that it does not apply at low momenta kj_ ^ Qs, which is the interesting 
regime for saturation. (6) In the second case, at k± <C Qs, our approximate solution 
is still a Gaussian, but with a different structure. It is scale-invariant and universal 



^Throughout, we consider the hadron in its infinite momentum frame, and use Ught-cone vector 
notations, v^^ = i;^, v^), with w+ = (1/V2)(t;° + v"^), = (1/V2)(u° — v^), and v_l = {v^,v^). In 
these notations, the hadron four-momentum reads = (P+,0,0_l), with P+ oo. 
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(up to logarithmic corrections which enter via the matching with the solution in the 
high-momentum regime), and has the right properties to describe gluon saturation. The 
associated gluon density increases only logarithmically with l/k'j_, a result formally similar 
to that obtained in the classical MV model [0, |l8l, but which here arises in a quite 
different way: The non-linear effects leading to saturation are encoded in the quantum 
evolution, and not in the solution to the classical equations of motion. Our results appear 
to be consistent with other recent analyses of the non-linear gluon evolution P, . 



In the MV model 0, gluon correlation functions at small x are obtained as classical 
correlations within the stochastic Yang-Mills theory with the following equation of motion 

(D^F'^'^Ux) = 6^+paix) , (1) 

where pa{x) is the effective color charge density at the "soft" scale k'^ = xP^ due to 
the "fast" partons with longitudinal momenta p'^ ^ k~^. Because of this separation of 
scales, the source is time-independent and localized near the light-cone, within a distance 
Ax~ ^ 1/k^. It is furthermore a random quantity whose spatial correlators can be 
summarized as a gauge- invariant weight function [p] , which is the probability for having 
a color charge distribution with density Pa{x) (with x = (x~,x_|_)). We have introduced 
here the momentum- space rapidity r = ln(P+/A;+) = ln(l/x) to indicate the dependence 
of the weight function upon the soft scale k^. 

Thus, equal-time gluon correlation functions at the scale k~^ = xP^ are obtained as: 
{A\{x^,x)J^,{x+,y)---)r = Jvp WAp]A'MX{y)-- - , (2) 

where A^^ = A\[p] is the solution to eq. (|1|) in the light-cone (LC) gauge A'^ = 0, which 
is the gauge which allows for the most direct contact with the gauge-invariant physical 



quantities |T^, |T^. For instance, the gluon distribution function is obtained as [|n 



xG{x,Q') = lj§^,^iQ'-kl){\K\k)\')^, (3) 

— * 

where J-"^^ = d~^Al^ is the classical electric field, and k = (A;"'",k_L) with k^ = xP^ = 
P+e""^. However, to have explicit expressions for these classical fields, it is preferable to 
express the LC-gauge solution A^J^p] in terms of color source p in the covariant gauge 
^^'A^ = (COV-gauge). One then obtains (with p = paT", etc.)lll: 



A'{x) = -U{x)d'U\x), 
9 

U\x~,Xi_) = Pexpi^ig J dz~ a{z~ ,x±)^ , 

-Vla(f) = p(x), (4) 

where pa, A^ = S^^Ua, and = —d^aa are respectively the color source, the vector 
potential, and the electric field in the COV-gauge. In the saturation regime, the source 
and the fields are strong, A^ ^ p ^ 1/g, and the classical problem is fully non-linear. 
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The source pa and the associated weight function Wr[p\ are constructed in perturbation 
theory, by integrating out the fast quantum gluons in layers of to leading order in 
asln(l/x), but to all orders in the strong background fields generated at the previous 
steps [l^, |15|, |l^. With the gauge-fixing prescriptions advocated in Refs. |jlO|, |l6l, this 



implies that the source Pa{x) has support only at positive x~, with Xq ^ x~ ^ l/fc"*", where 
Xq = 1/ is the Lorentz-contracted longitudinal size of the hadron. By introducing the 
space-time rapidity y = ln(a;~/xg) and recalling that = l/{xP^) = Xgc"^, we 

conclude that pa has support at < y < r. 

The evolution of the weight function WV[p] with increasing r = ln(l/x) is described 
by a functional renormalization group equation originally proposed by Jalilian-Marian, 
Kovner, Leonidov and Weigert and explicitely constructed in Refs. |l^ (see also 



Refs. p2|). It acquires its simplest form when written as an equation for PVt-[q^] (recall 



that a and p are linearly related, cf. eq. (^), in which case it reads flO 

dWJa] 



dr 

with the following, positive definite, Hamiltonian: 



-HWr. (5) 



J 2n {z_i_ — x_i)^ oa'^{x±) 
This is the same Hamiltonian which generates Balitsky's evolution equations, as rewritten 



in functional form by Weigert fTj] ; this establishes the equivalence between the approaches 
of Refs. § and |TD|, |T§. 

In the equations above, 

= V\x±_) = U^ix" =Xoe^,x_L), a°(x^) = a"(x" =x^e'^,x±), (7) 

are the Wilson line and the COV-gauge "Coulomb" field in eq. (^ evaluated at x~ = 
XqC"^ = l/k"^, which is the largest longitudinal coordinate for the color source at the 
scale = P~^e~^ . Thus, the functional derivatives in eqs. (0)-(P) are to be taken with 
respect to the color field in the highest bin of (space-time) rapidity t < y < r + dr, where 
the quantum corrections are located. This suggests a space-time picture of the quantum 
evolution where the momentum rapidity r and the space-time rapidity y are interrelated: 
The classical source 

Py{x±) = P"(X" = XqB^, x±) (8) 

is constructed in layers of space-time rapidity, with the contribution in the rapidity bin 
(y, y + dj) obtained by integrating out the quantum gluons with longitudinal momenta 
in the momentum-rapidity bin (r, r + dr) with r = y. This relation will become explicit 
in the solutions to the RGE to be found below. 
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Eq. d^) is a functional Fokker-Planck equation with "time" r. It depicts the quantum 
evolution towards small x as the diffusion of the probability density Wr [a] in the functional 
space spanned by (ya{x~ ,x±). Since the r.h.s. of this equation can be written as a total 
derivative with respect to a [0, |T^, the correct normalization: 

Va Wr[a] = 1, (9) 

is automatically preserved by this evolution. In the weak field limit, where the Wilson 
lines can be expanded to lowest order in a, eqs. ®-® reproduce the BFKL equation for 
the charge-charge correlator {p{k±)p{—k±)) |T3|, In the strong field regime a ~ 1/g 



specific to saturation, eq. @ is equivalent to an infinite hierarchy of coupled evolution 
equations for the correlators of a. In what follows, we shall not study these coupled 
equations, but rather attempt to solve directly the functional RGE (|]) in some limiting 
cases. More details will be presented somewhere else . 



As announced in the introduction, we shall find approximate solutions to eq. (j^) in two 
different kinematical regimes: k±_ ^ Qs{'t) and <^ Qs{t), where k± is the transverse 
momentum at which we measure the correlation functions (i.e., the transverse resolution of 
the external probe), while Qs{t) ("the saturation momentum") is an intrinsic momentum 
scale which is introduced by the initial conditions (note that there is no explicit scale in 
the evolution equation (§)-(|^)) and which increases with r. No attempt will be made to 
describe the intermediate behaviour at k± ~ Qs, and thus the onset of saturation. Because 
of that, the saturation scale itself will be not exactly determined, but only estimated via 
a study of the onset of non-linearity in the high momentum regime. 

a) The high momentum regime. As we shall shortly discover, at least within 
the approximations to be considered here, it is the same scale Qg which controls both the 
transverse correlation length in the problem, and the onset of the non-linear regime. That 
is, the typical variation scale for the Wilson line V{x±) is 1/Qs, and the non-linear effects 
become important at transverse momenta of order Qg or less. Thus, as long as we are 
probing the system on a much shorter resolution scale, corresponding to some external 
momentum k± ^ Qs, we can neglect the non-linear effects and perform a short-distance 
expansion in the Hamiltonian (|^) : 

1 - V^V^ ^ {z'J - x^){&V^)^^^ ^ ig{z^ - x^)d^a{x±), (10) 

where (with x~ = x^e'^) 

a''{x±) = [ dx~a°'{x~ ,x±) = [ dja^{x±), (11) 







is the effective Coulomb field in the transverse plane, and we have introduced, for further 
convenience, the notation (x_l) = x~ a'^ {x~ , x j_) for x~ = XgC^. 



With these approximations, and after simple manipulations to be detailed in Ref. 
the Hamiltonian (P) is brought into the following form: 
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The RGE associated to this Hamiltonian is still non-linear, and we have not been able to 
solve it exactly. To make progress, we perform a mean field approximation in which we 
replace, within eq. (p!^), 

d'a''\x±)d'a'"'{yi_) {d'a''''{xi_) d'a^%y±))r = N6^^V%{x± - y±), (13) 

where the expectation value in the r.h.s. is defined as in eq. (0), and we have assumed 
homogeneity in the transverse plane, for simplicity. (The color structure in the r.h.s. 
follows from gauge invariance.) With this approximation, the RGE is linear and can be 
solved by a Gaussian weight function. The correlation function in eq. ([T3|) can be 
then computed in terms of the width of this Gaussian, which, as we shall see, entails a 
self-consistency condition describing the evolution of the width with r. 

Specifically, we have to solve the following equation: 



^ J d^x± J Sy±_ Vr{x± - y±_ 



dr 2 J -V - ""^^ 6a^{x±)6a^{y^)' 

g^N 1 

'Drix±-y±) = {x±\=^\y±)^r{x±-y±), (14) 

vr V_L 

or, in momentum space. 



where in writing the approximate equality we have used the fact that k_i_ is the hard 
external momentum, while p± is relatively soft, p± ^ Qs kj_. Since fc^a,-(fc±) = Prik^), 
cf. eq. (0), the Hamiltonian in eq. ([l^) is finally rewritten as (with as = g'^/Air): 

n ^ ( ^'^^ r (1. 



2 J (27r)2-"^ ' Sp-ik^)6p-{-ks_y 

k r}^ 

G^k^) ^ Aa.N r -^pIUpl). (16) 



(27r)^ 

Thus, in this high-momentum regime, the solution to the RGE is most naturally written 
as a Gaussian in pa{x). The only subtle point about this solution refers to its longitudinal 
structure. To understand this structure, note that the two functional derivatives in eq. ([T6| ) 
act at the same point x~, namely at X — Xj. — Xq e^. This implies that the correlations 
generated by the RGE are local in x~ . We thus search for a solution PVV[p] of the formQ: 

"2/0 -^2py^k^)K\YM)p;{-k^)Y (17) 

where we have used the space-time rapidity y to indicate the ^"-dependence of the various 
functions and defined Py{x_C) = x~p''(x~,x^) for X X Q e^. The integration over y in 

^The initial condition for the solution will be discussed in Ref. 1231. 
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eq. ( p!7D runs effectively up to r, but the upper integration limit is conveniently absorbed 
in the support of the kernel A~^(y, kj_) (see eq. (^) below). The overall factor AC follows 
from the normalization condition (^) as (up to some irrelevant, r- independent, factor): 

K = [det A.(y, k^r'/' = n n [^^(y, k^T^'''~'^^\ (18) 

kj_ y 

where in writing the second equality we have considered a lattice version of the 3- 
dimensional configuration space, with discrete points k± and y. 

It is now straightforward to verify that the functional (|l^) satisfies the RGE associated 
to -f^high-fex provided the width A^(y, k_i) obeys the following equation: 

^^l|M = i(,_y)(;,(ij, (19) 

which shows that the evolution of the width with the momentum- space rapidity r takes 
place at the space-time rapidity y = r, that is, at the end point of the color charge 
distribution at the previous step. The solution to eq. (|19|) is immediate: 

A.(y,A;^) = 9{T-j)g^{k^). (20) 

This gives the following 2-point function for the color charge: 

{P;i^±) Py'iyi.))r = 6^'5iy - y')e{r - y)g,{x^ - y^). (21) 

Note that, for y < r, the correlation function Gy{x± — y±) is independent of r, since 
determined uniquely by the quantum evolution up to the momentum rapidity y (via 
eq. (|25|) below). 

According to eqs. (0) and (0), the width of the Gaussian is the inverse of a 6'-function, 
which means that the weight function Wr[p] is identically zero for all the functions Py{x±) 
having support at rapidities y > r. Thus, in the functional integral (^ one can freely 
integrate over all the functions Py{x±), without any restriction on their support (other 
than y > 0); the restriction to y < r will be automatically taken care of by the weight 
function. 

At this point, one should recall that the correlation function eq. (0), and thus 
the function Qrik^) in eq. (0), are expectation values with the weight function ([T^), and 
thus functionals of A,-. Thus, eqs. (|19D or (pOD are really self-consistency conditions, to be 
made explicit now. To this aim, it is convenient to integrate eq. (^) over y, 

= Qr{k^), (22) 

and thus obtain an evolution equation for the quantity 

Pr{kl) = I djXr{j,k^), (23) 
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which, in the present weak field regime, is the unintegrated gluon distribution. To see 
this, note that, to lowest order in p, JF+J' ~ i{k^ /k\)pa , so eq. (H) becomes 



TT 



(27r)2 kl ^1 /r 47r 7o kl 



2 



where is the hadron radius, and we have also used eqs. ( ^TD and (|23|) . 

From eqs. ( [1^ ) and (|2T|)-(p3D, one obtains ^r(p±) = i)- / vV) I^t{j>i) i which, together 
with eqs. (|1^) and (p2D, provides an evolution equation for Urik^) : 

^^^^^"^^ - /..(Pi) . (25) 



dT TT 7o 



By also using eq. (p^), this can be recognized as the evolution equation for the gluon 
distribution function in the double logarithmic approximation (DLA) : 



■xG(x,g') = ^xG(x,g2). (26) 



Given the approximations that we have performed, this is indeed the expected limit of 
the evolution equation [^. The solution to eq. (|26|) is well known ||2^ : if one holds as 
fixed (independent of Q^), the solution grows like: 



xG{x, Q'') oc exp <( 2 r HQyQl) } . (27) 



For our purposes here, what really matters is that the initial condition to eq. ( pSf ) intro- 
duces a momentum scale in the problem, while there was no such a scale in the RGE 
by itself. For instance, in a simple valence-quark model, which should be a reasonable 
approximation near a; ~ 1 (or r ~ 0), one has [|19| (with Cp = [N"^ — 1)/2N) 



2 



xG{x, Q2) - '^lH^L In for a; ~ 1 , (28) 

^QCD 

where ^'qcd enters naturally as an infrared cutoff when taking into account the over all 
color neutrality of the hadron on scale sizes of order 1/ Kqcd [H|- Thus, for r ~ 0, 



/io(fci) ~ 2as/R\ (29) 

which for a proton is a relatively small momentum scale to start with, but which at low 
X is enhanced by the quantum evolution described by eq. (p5|). 

When r, and therefore /Xt, are sufficiently large, the fiuctuations of the color field 
described by eq. (0) become large as well, and the non-linear effects start to play a role. 
To study the onset of non-linearity, and verify the approximations performed in eq. (p!0D , 
it is useful to consider the 2-point function {V^Vy)r- This is easily computed by expanding 
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the path-ordered exponentials within the Wilson lines, performing the contractions of the 



fields a with the help of (cf. eqs. (0) and (p2D ) 

1 <9^r(P±) 



pi 



(30) 
(31) 



dr dr 

and then recognizing the result as the expansion of an ordinary exponential. (See Refs. 



T7| , |T^ for similar calculations.) The final result can be written as: 



{V\x^)V{0^))r = exp -g^N / dy 7^(0^) - 7y(xx) 



with 



d'^P± fJ,r{p] 



P\ 



Since /ir(pi) is rather slowly varying as a function of p±, the integral in eq. 
dominated by small momenta ^ V^l) where we can approximate 



f d'^p_L 1 






/ (27r)2 p\ 







rfV ^ {P± ■ a^±)' 
(2vr)2 pi 



xl 
~ m 



2 IGvr 

By performing a similar approximation on eq. (BSf), one obtains: 



r2 A2 



(32) 

(33) 
i is 

(34) 



(^nxx)^(Ox)), 



exp 



xl 



PI 



2 /^r(pi) 



(35) 



By inspection of this equation, it should be clear that, as anticipated, it is the same scale 
Ql ^ l/xl which controls both the non-locality and the onset on non-linearity: this is the 
scale for which the exponent in eq. ( P^D is of order 1. This implies the following estimate 
(actually, an equation) for the saturation scale: 



(36) 



For k± ^ Qs, the approximations in eq. (|T0D are justified, and the Gaussian (|T7D is the 
correct solution to the RGE within the mean field approximation. But this solution 
cannot be extended at low momenta k± ^ Qg. This shows the inconsistency of some 
previous calculations within the classical MV model, where the Gaussian weight function 
(|1^) has been used at all momenta, including in the saturation regime [|1^, |1^ . 

For later reference, it is useful to combine eqs. (p6l), (p4|) and (|27| ) to obtain an estimate 
for the r-dependence of Qs in the DLA (with a<j = asN/ir) : 



QI{t) oc xG{x,Ql) ~ exp 2Va.rln(Q2/Q2 



or QI{t) oc e 



4o.r 



(37) 
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b) The low momentum regime. We now turn to the more interesting regime at 
small momenta, k_i_ -C Qs{t), when the hadron is probed over distances large compared 
to the correlation length 1/Qs{t) (but still small as compared to 1/Aqcd)- In addition 
to being relatively rapidly varying, the electric fields in this regime are also expected to 
have large amplitudes, of order 1/g (see eq. ( |5T| ) below). In this regime, operators like 
V^Vx are rapidly averaging to zero. Thus, in a first approximation, we shall simply ignore 
all the Wilson lines in the Hamiltonian (^). We shall verify a posteriori that this is a 
consistent approximation, by computing the expectation value {V^Vy)r with the resulting 
weight function. This approximation will be further justified in p3[| . 



In this "random phase approximation" (RPA), the Hamiltonian reads simply: 



and does not involve the coupling constant at all. Thus, this is formally like a free theory, 
although it has been obtained in a strong field regime where ga ^ 1. In fact, this is 
not really a free theory, since the classical fields whose correlators are needed are highly 
non- linear functionals of ga, cf. eq. @. However, as we shall shortly see, these classical 
non-linear effects are not essential for getting saturation. 

To determine the weight function one has to integrate the RGE (|^) over all 

rapidities r' with < r' < r. Clearly, the RPA leading to the Hamiltonian (|38|) , which 
requires kj_ -C Qs{t'), is not appropriate for all such intermediate r'. Let f{kj_) be the 
rapidity at which the saturation momentum becomes equal to the external momentum^ : 

Qlifik±)) = kl. (39) 

Then eq. ( ^8]) applies at r' > t(A;^), while for < r' < f{k±) we are rather in the 
high-momentum regime. This suggests the following approximation: 

Wr[a] ^ W^^\a]W}^ZrW\. (40) 

where IV^ '^'^[a] is the weight function in the high- momentum regime, eq. (|T7|), evaluated 
at f{k_[_) (this is the initial condition for the evolution described by ifiow-fcx)) ^^nd 
is obtained by integrating the RGE with Hamiltonian ( ^8]) from f{k_[_) up to r. 

A crude estimate for the separation scale f{ki_) can be obtained by using the solution 
in the high- momentum regime. Eq. ( pTf ) implies Q1{t) oc e^"*"^, where c = 4 in the DLA 
but will be kept here as a free parameter, to partially account for our ignorance of the 
true dynamics at /c_l ~ Qs- This, together with eq. (|39D, implies 



^We are grateful to Al Mueller for pointing us the importance of the separation scale t(A;j 
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which shows that, when k± ^ Qs{t), we have also r ^ f{k±). Physically interesting 
correlation functions like the gluon distribution (|^) involve generally an integral over all 
rapidities < y < r. In the regime r ^ ^^(^x), the dominant contributions to such 
correlations come from the interval f{k±) < j < t, and are thus determined by the piece 
of the weight function (|40D. Given the simplicity of the Hamiltonian (pSj), this 



is easily obtained as |^ 

= AT.exp |- i Q^^ dy J d'x^ d'a^{x^) C\y) d'a^ix^)"^ , (42) 
with the following width: 

C.(y) = {l/n)e{T-y). (43) 
This implies (for f{k±) <j,j'< r) 

{a^{x^)a;,{y^))r = {lM6^'6{j - j') {x^\^\y^), (44) 
{d'a^{x^)d'a^,{y^))^ = {l/n)6^''6{j - j') 6^'\x^ - y^), (45) 



showing that the probability distribution for the COV-gauge electric field = —d^a is 
local and homogeneous both in the transverse space and in space-time rapidity (within 
the interval f{k±) < y < t). This is consistent with our assumption that transverse 
correlations in the system are over a typical scale 1/Qs, and thus cannot be resolved 
when the hadron is probed with a much lower resolution. But this also shows that the 
transverse 5-function in eq. (^) must be taken with a grain of salt: this is really delocalized 
over the correlation length 1/Qs{j). 

It is now possible to compute the gluon distribution function in this low momentum 
regime. Since the fields are strong, we have to use the fully non-linear expression for the 
classical field in eq. (^). Eq. (^ involves the following 2-point function: 

(^+^(fc+,x^)^+^(-fc+,l/^)), = J dx- j dy-e^'^^^--y~\:Ft\x-,x^):Fi\y-,y^))r, (46) 

which we compute as follows: 

= {d^a\x)d^a\y))^{ul{x)U,M)^, (47) 

where, as indicated in the second line, the two COV-gauge electric fields d^a can be 
contracted only together (mainly because of the path ordering in the Wilson lines, which 
forbids other contractions [|17[ ^). By also using eq. (^Sf ), one obtains: 

{:Fr{k\x^):F^\~k\y^))r ^ / dy5^'\x^~y^) {ul{x^)U,{y^)) . (48) 

TT (fcx) 
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The transverse (5-function in the r.h.s. of eq. ( P8| ) imposes x± = y_L, and therefore 
{UjUy) 1. More precisely, this 5-function is spread over a distance \x_i — y±\ ~ 1/Qs(j), 
which is precisely the correlation length for the 2-point function {U^{x±)Uy{yj_)) . Over 
such a distance, the latter decreases by a factor b, with 6 > 1 but not much larger (typi- 
cally, 6 ~ e). We thus obtain: 

{J'rik\x^)J^^%-k\y^))^ = ^:^l-lls(^)(x^-y^){r-fik^)). (49) 

Up to the global factor 1/6, this is the same result that would have been obtained in the 
linearized, or weak field, approximation JF+* ~ —d^aa- AH the non-linear effects — which 
were a priori encoded in the 2-point function of the Wilson lines — have dropped out 
from the gluon distribution because of the locality of the propagator (^). 

By using eqs. (^ and (^91) , one can compute the number of gluons per unit of transverse 
phase space in this low-momentum regime. One obtains {b± is the impact parameter) : 



(Pk I (Pb I 



/ -£^^"'^"'H^r{k\x^)Tr{-k\Q^))r = -^^i{^-^ik^))^ (50) 



or, after also using eq. (|T|) 



<P(xG) _iV^-l 1 i^QM. (51) 



where the two unknown constants b and c have been combined into a = be. 

Eq. (^) is one of the main results in this paper. It shows marginal saturation (when 
k'j_ decreases, the gluon density still increases, but only logarithmically), and is consistent 
with the unitarity bounds: at fixed k^, the gluon density increases only linearly with r 
(this is most obvious on eq. (^0])), and so does also the associated distribution function 
xG{x, Q^) for Q'^ <ti Q^. The latter is easily obtained from eq. (pl|) as (with / (Pb± = nR"^): 

n"^ — ^ ^ r 

xG{x,Q') = ^^_i?2g2 ln(Q^(r)/g') + l 

/V^ — 1 1 r 1 

= — -rp-W ca.(r - f(g)) + l . (52) 

If extrapolated up to Q ~ <5s, this agrees quite well with the corresponding extrapolation 
of the gluon distribution at high momenta, as given by eqs. ( plD and (|36|) . 

According to eq. (|5l|), the gluon density at saturation is of order l/a^, which corre- 
sponds to color fields as strong as ~ 1/(7. Still, as the above analysis clearly shows, 
the non-linear effects leading to saturation are not those in the classical field, but rather 
those in the quantum evolution, which led to the local weight function (|4^). In other 
terms, the saturation is built-in in the effective action at low momenta. In this respect. 



our conclusions differ from those in Refs. |17, 18], although the final results found there 



for the gluon density are formally similar to our eq. (^ 
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On the other hand, our results here are consistent with some recent analyses of the non- 
linear gluon evolution by Mueller and Levin and Tuchin [jl2[. In Ref. [jl2|, Levin and 



Tuchin have obtained exact, semi- analytic, solutions to the Balitsky-Kovchegov equation 
for {V^ {x±)V{y±))r P, H, and then used these solutions to estimate the gluon distribution 
function. Their result is consistent with our eq. ( p2D provided one takes c = 4, which 
coincides with the prediction of the (a priori, crude) DLA (cf. eq. (pT])). 

As a check of the self-consistency of the RPA, let us compute the 2-point function 
{V''{x±)V{Oj_))r for x]_ = l/k\ > 1/Q1{t). This is again given by eq. ( p2D where the 
integral over y is now decomposed into two pieces: < y < f {%_[_) and f{xi_) < y < r, 
with f{x±) = f{k'j_ = l/x]_) (cf. eq. (|39D). The first rapidity interval gives an attenuation 
factor b e, while in the second interval we can use the propagator jyixx — y±) from 
eq. (|4^). This gives: 

7.(0x) - Jy{x^) = ^ / 1^ ^ [l - ^'^-^^ - ln(g^(y)4) , (53) 

where the scale Qs(y) has been introduced to cut off an ultraviolet divergence. It can be 
shown [pi], via an analysis of the Balitsky-Kovchegov equation, that this is indeed the 



natural UV cutoff introduced by non-linear effects. Since (cf. eq. (^if )): 

HQl{y)xl) = m,(y-f(x^)), (54) 

we finally obtain: 

{V\x^)V{0^))^ ~ iexpj-ca^l^^^ rfy (y - f (xx))| (x exp{-(c/2)«^(r - f (x^))^} 



oc exp 



It{xj_) 

^ HQl{T)xl) - \ . (55) 



2c 



Once again, this coincides with the corresponding result in |T2[ provided one takes c = 
4. Eq. ( ^51) shows that the correlator of the Wilson lines is rapidly decreasing when 
QK''")^! ^ so that the RPA is indeed justified, at least as a mean field approximation. 



Note finally that, in contrast to the weight function at high momenta, eq. {\17\), the 
width of the Gaussian in eq. (^2|) does not involve any mass scale. Thus, the weight 
function at low momenta is scale invariant and universal (in the sense of being insensitive 
to the initial conditions, and therefore the same for all hadrons). These properties trans- 
mit to the correlation functions computed with this weight function, up to logarithmic 
corrections which enter via the lower rapidity limit f{k±) (which depends logarithmically 
on kjj and via the ultraviolet cutoff Qs{t). Thus, the correlation functions at saturation 



are functions of hi{Qg{T)/k\), as manifest on eqs. (pi]) and (|55D . 



To conclude, in this paper we have constructed approximate solutions to the RGE for 
the Color Glass Condensate valid at either large, or small, k^ momenta. In both limits, 
we have used a mean field approximation. There is no difficulty of principle for using a 
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mean field approximation in more general cases, since this inevitably leads to a Gaussian 
form for the functional For strong fields, the width of this Gaussian will satisfy a 

non-linear evolution equation, which in the high-momentum regime can be shown to 
be consistent with a previous result by Mueller and Qiu 0. Beyond the mean field, we 
would like to be able to understand the limitations of the mean field approximations and 
perhaps to fully solve the renormalization group equation. 
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